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ABSTRACT 

The Sunyaev-Zeldovich (SZ) effect is a promising tool to study physical properties of 
the hot X-ray emitting intracluster medium (ICM) in galaxy clusters. To date, most SZ 
observations have been interpreted in combination with X-ray follow-up measurements 
in order to determine the ICM temperature and estimate the cluster mass. Future 
high-resolution, multifrequency SZ observations promise to enable detailed studies of 
the ICM structures, by measuring the ICM temperature through the temperature- 
dependent relativistic corrections. In this work we develop a non-parametric method 
to derive three-dimensional physical quantities, including temperature, pressure, total 
mass, and peculiar velocities, of galaxy clusters from SZ observations alone. We test 
the performance of this method using hydrodynamical simulations of galaxy clusters, 
in order to assess systematic uncertainties in the reconstructed physical parameters. 
In particular, we analyze mock Cerro Chajnantor Atacama Telescope (CCAT) SZ 
observations, taking into account various sources of systematic uncertainties associated 
with instrumental effects and astrophysical foregrounds. We show that our method 
enables accurate reconstruction of the three-dimensional ICM profiles, while retaining 
full information about the gas distribution. We discuss the application of this technique 
for ongoing and future multifrequency SZ observations. 

Key words: cosmology: observations - galaxies: clusters: general - X-rays: galaxies: 
clusters - cosmic microwave background 



1 INTRODUCTION 

Clusters of galaxies are the largest gravitationally bound 
structures in the universe, promising to serve as powerful 
cosmological probes. In the current hierarchical structure 
formation paradigm, galaxy clusters form out of the grav- 
itational collapse of large peaks in the primordial density 
fluctuations in the early universe and grow through mass ac- 
cretion. Since their evolution traces the growth of the linear 
density perturbations, galaxy clusters are powerful laborato- 
ries f or studying the na t ure of dark matter an d dark energy 
(e.g., lAllen et al.l |200S| ; IVikhlinin et al. 1 120091 ; iMantz et alj 
I2010I ). 

Galaxy clusters can be detected through X-ray emis- 
sions from the hot (0.5 < T < 10 keV) intracluster medium 
(ICM). In recent years there has been significant advances 
in discoveri es of galaxy clusters via th e Sunyaev-Zeldovich 
(SZ) effect jSunvaev fc Zeldovich|[T97oh - a spectral distor- 
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tion of the cosmic microwave background (CMB) due to 
the inverse Compton scattering of the CMB photons off 
by the hot X-ray emitting plasma in clusters. The SZ ef- 
fect signal is independent of redshift, making it an excel- 
lent method for finding clus ters in the high redshift uni- 
verse l|Carlstrom et al] 12002). Moreover, since the SZ in- 
tensity depends linearly on the density, unlike the X-ray 
flux, which depends on the squared density, the SZ sig- 
nal is inherently sensitive to the global properties of clus- 
ters out to and beyond the virtualized regions of clus - 
ters (|Plagge et alj|20ld : IPlanck Collaboration et~ai1 |2012|) . 
Dedicated multifrequency SZ observations, such as ACTq 
Planck 0, and SP Tq are discovering a large sample of high- 
redshift clusters (IVanderlinde et al.l l2010t iMarriage et al.l 
l201ll : IPlanck Collaboration et al.l l201lh .~ promising to pro- 
vide new insights into the cluster gas physics and cosmology. 



http: / /www. phy.princeton.edu/act 
http://www.rssd.esa.int/index.php?projcct=planck 
http: / / pole.uchicago.edu 
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At the same time, a number of high-resolution SZ exper- 
iments including ALMifl CARMAQ CCAT0, and MUS- 
TANQj, are also underway or planned, promising a dramatic 
increase in sensitivities, spatial resolution, and spectral cov- 
erage. Current high-resolution SZ observations of individual 
clusters have already revealed rich phenomena (e.g., shocks, 
substructures, relativistic particles) i n the hot gaseous at- 



mospheres of merging clusters (e.g., Komatsu et al. _ 200 ll; 

^' 20ld;IColafrancesco et alj|201ll;lKorngut et all 



Maso n et al 
20 111; 



Menanteau et al 



2012 



Mroczkowski et alj l20ll 



Zemcov et al.ll2012i ; IProkhorov fc Colafrancescdl2012f ) 



One of the fundamental limitations of the current SZ 
effect observations lies in their inability to measure the ICM 
temperature, independently of the X-ray measurements. 
Knowledge of the clusters temperatures is especially critical 
for estimating the cluster mass by assuming the hydrostatic 
equilibrium of gas in the gravitational potential. Cosmolog- 
ical constraints from the SZ surveys thus require extensive 



ray satellites (e.g., Reese et al. 2002 


; LaRooue et alj|2006l: 


iBonamente et al.ll20od; Benson et al 


1201 lh. 



Over the next years, high-resolution, high-sensitivity SZ 
observations will enable detailed studies of the thermody- 
namics structure of the ICM and mass modeling using the 
SZ data alone. The idea i s to use the temperature-dependent 
relativistic corrections dWrighd 1 19791; Is azonov fc Sunvaevl 
1 19981 ; lltoh et al.|[l99g| ; IChluba et alj|2012h that are signifi- 
cant for a hot (Tx J> 3 keV) galaxy clusters. The relativistic 
corrections introduce additional frequency-dependent sig- 
nals, making it possible to derive the gas temperature from 
mult i frequency SZ observations, independently of X-ray 
data llPointecouteau et~a l ll99j; lHansen et alJEool : iHansenl 
12004 IProkhorov et al] 1201 ll ). Upcoming SZ experiments 
should further enable measurements of the peculiar veloc- 
ity of the cluster and internal bu lk, turbulent, and rotational 
gas motions within clusters (e.g.JChluba fe M annheim 2002; 
iNagai et al.1l2003l ; ISunvaev et al.ll2003l ; iDiego et alj|2003h . 

A variety of techniques have been developed to derive 
three-dimensional (3D) profile of physical quantities from 
two-dimensional physical observables on the plane of the 
sky. These techniques can be divided into two broad classes 
that will be denoted here as inverse method and forward- 
fitting method. The forward-fitting approach assumes a pri- 
ori parameterized models for the 3D physical quantity in 
order to compare it s projection with the observables in 
the p la ne of the sky jMroczkowski et al]|2009l ; iBulbul et al] 
l2010h . IColafrancesco fc Marchegianil (^PloT T for example, 
have developed a formalism to derive the 3D ICM tem- 
perature profile from spectroscopic SZ observations using 
a parametric model. However, one of the downsides of the 
forward-fitting method is a need for an a priori parame- 
terization, which might not wholly reflect the 3D physical 
quantities and not adequately retain the full information 
e.g. of the temperature profile in the gas distribution. Fur- 
thermore, the desired physical parameters of clusters ex- 
tracted in this way are prone to large systematic uncer- 
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tainties, due to simplified modeling of the properties of 
the hot ICM and cluster morphology. The inverse method, 
on the other hand, introduces one or more arbitrary pa- 
rameters to guarantee a smooth and physically acceptable 
3D physic al quantity by combinin g X-ray and SZ imag - 
ing data l|Yoshikawa fc Sutol Il999l; lAmeglio et all 120071) . 
Using a parametric model, IColafrancesco fc Marchegianil 
(|201dh have further developed a formalism to derive the 3D 
ICM temperature profile from spectroscopic SZ observations 
alone. 

The main goal of our work is therefore to develop a 
robust non-parametric approach to derive 3D physical pa- 
rameters of galaxy clusters from multifrequency SZ obser- 
vations. In order to assess how accurately one can recover 
the 3D ICM profiles, we analyze mock CCAT observations 
of simulated galaxy clusters extracted from hydrodynam- 
ical simulations, taking into account various instrumental 
noise and astrophysical contaminations. We also briefly dis- 
cuss the application of our technique for upcoming multi- 
frequency SZ observations. Throughout this work we as- 
sume the flat ACDM model, with matter density param- 
eter Q. m = 0.3, cosmological constant density parameter 
Qa = 0.7, and Hubble constant H = 100/ikms^ 1 Mpc -1 
where h — 0.7. Unless otherwise stated, quoted errors are at 
the 68.3% confidence level. 



2 THE SUNYAEV-ZELDOVICH EFFECT 

The SZ effect is a small distortion of the spectra of the CMB, 
due to the inverse Compton scattering of the CMB photons 
off of hot (~ 10 7 - 10 s K) electrons of the ICM trapped 
into the gravitational potent i al well of the dark matter halo 
l|Sunvaev fc Zeldovichl [l970l ; iBirkinshawl Il999h . The CMB 
photons have a low probability (r ~ 0.01) of interaction 
with the hot electrons of the ICM, which increases on aver- 
age their energy approximately by a factor O = fceT/meC 2 , 
producing a distortion of the black-body spectrum of the 
CMB: this appears as a decrement of the monochromatic 
flux of the CMB at frequencies smaller than 218 GHz and 
as an increment at frequencies larger than 218 GHz. 

The SZ effect is expressed as a small variation in the 
brightness intensity AI(v) of the CMB as a function of the 
observation frequency: 

AIM 



(It 



la 



• J P(r)*(i/;T(r),0(r))d* 



(1) 



where ar is the Thomson cross-section, P(r) = n e (r)fcf,T(r) 
is the electron pressure of the ICM at the volume element 
of coordinate r, kt is the Boltzmann constant, z is the line 
of the sight distance, I = 2(fc b T cmb ) 3 /(/ic) 2 , T cmb = 2.725 
K, and f3 — v p /c, where v p is the streaming velocity of the 
IC gas along the line of sight (positive if the IC gas is ap- 
proaching the observer). 

^{v; T(r), /3(r)) takes into account the spectral shape 
of the SZ effect and it reads: 



*(v;T(r),0(r))=g(v;T(T)) 



O(r) 



fc(i/;T(r),/9(r)) . (2) 



The first term in Equation @ is the thermal SZ (tSZ) dis- 
tortion with spectral shape: 

f + 1 



giy\T) = 



(e* - l) 2 



1 



4 (1 + o g (x;T)) , (3) 
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where x = hv/kT CTO h accounts for the frequency dependence 
of the SZ effect, an d for the relativis tic corrections related to 
the term o g (x,T) ()ltoh et al.lll998T h Its magnitude is given 
by the Comptonization parameter: 

y=^[ P ( r )dz, (4) 

m e c^ J 

If we assume a simplified isothermal parameterization of the 
ICM we have y = TkT/(mc 2 ), where r = <tt / n e dz is the 
optical depth. 

Note that the relativistic correction term o g (x,T) in 
Equation (O becomes significant at high frequency and/or 
high ICM temperature: for a hot Tx > 7 keV clusters, the 
relativistic corrections are ~ 5% of the total signal at 100 
GHz and ~ 30% at 675 GHz, making the SZ brightness 
very prone to the ICM temperature. The second term in 
Equation ([2]) is the kinetic SZ (kSZ) distortion, which arises 
from the Doppler shift due to bulk motion of the electrons 
with respect to the rest frame of the CMB photons. The 
spectral shape of the kSZ effect is given by, 

h(u;T,/3) = { ff 1)2 (l + o h (x;T,f3)) , (5) 

where the ter m Oh(x;T, 0) is th e relativistic corrections of 
the kSZ effect (|ltoh et al.lll998T l. and its magnitude is pro- 
portional to ;8/0. 

It is worth noting that in Equation ((2| the kinetic 
distortion is maximal at the crossover frequency (i>q = 
218 GHz), where the tSZ effect vanishes, while at lower or 
higher frequencies it is smaller than the thermal distortion 
of a factor ~ /3/0 ~ 10~ 2 — 10 _1 for typical values of gas 
temperature and peculiar velocity relevant for clusters. Note 
also that, in order to calculate g(u;T), an a priori measure- 
ment of the gas temperature e.g., from X-ray is needed. X- 
ray temperature is also essential for estimating the cluster 
mass through the hydrostatic equilibrium assumption. For 
these reasons, SZ studies of cluster have have been mostly 
done in combination with X-ray measurements to date (e.g. 
iLaRoaue et~alll2006h . 

We observe that, given the spectral signature of 
^(v;T(r), /3(r)), in principle it is possible to separate the 
physical variables P(r), T(r), /3(r) by deprojecting AI(y) in 
Equation Q, in order to infer P(r) *(i/;T(r),#(r)). Sen- 
sitive multifrequency SZ observations can measure AI(u), 
thereby enabling measurements of the desired physical pa- 
rameters, including P(r), T(r), /3(r), independently of the X- 
ray measurements. 

Specifically, this reconstruction is accomplished by com- 
paring the theoretical prediction of P(r) ^(v; T(v), j3(r)) to 
the deprojected observed quantities. In the next section we 
present a robust non-parametric approach to reconstruct 3D 
physical quantities from two-dimensional projected SZ maps 
by inverting Equation Q. 



3 SPECTRAL DEPROJECTION TECHNIQUE 

A variety of techniques have been developed to derive 3D 
physical quantities (e.g. P(r) ^(f, T(r), /3(r))) from two- 
dimensional physical observables on the plane of the sky 
(e.g. AI(y)). As discussed in §[T] in this work, we develop a 
robust non-parametric approach for inverting Equation {TJ 



and constraining 3D temperature and pressure profiles from 
SZ-only multifrequency observations. 

We assume that the cluster is spherically symmetric 
and it has a onion-like structure with n concentric spher- 
ical shells, each characterized by uniform gas density and 
temperature within it. Therefore, the cluster image in pro- 
jection is divided into m rings (or annuli) of area A = 
(Ai, A2, A m ), which are assumed to have the same radii 
of the 3D spherical shells. The SZ boundary (the most ex- 
ternal ring) will be represented by the mth ring of area A m 
and radius R m . Let us define £j as the SZ signal to be re- 
covered from the deprojection method within the ith shell. 
In our analysis s and e read: 

s = AJ(i/)/J„, e =-^P(r)*(i/;T(r),0(r)). (6) 

with r = Pi, i = l, m. In this way, the contribution of 
the i-th shell to the SZ surface brightness in the ring j of 
the image will be given by 

s = (V#e)/A (7) 

where the operator # indicates the matrix product (rows 
by columns). The matrix V- is a rectangular m x n matrix, 
with the m-dimensional column vectors V 1 , V 2 , V n rep- 
resenting the effective volumes, i.e. the volume of the j-th 
shell contained inside the i-th annulus (with j i) and 
corrected b y the gradient of e in side the j-th shell (see Ap- 
pendix B of lMorandi et al1l2007T ). We assumed that n > m: 
this means that in jth ring of the image we account also 
for the so-called edge effect, i.e. for the contribution of the 
shells [m + 1, m + 1, .., n) of radius R > Pbound- 

Now, in order to infer e, we introduce two types of in- 
formation. The first is a statistical information and allows 
us to obtain a solution which reproduces the data with a 
known accuracy. In our case, since the error on the data are 
assumed to follow a normal distribution, the x 2 gives the 
accuracy of the solution: 

X 2 = [s - (V#e) /A^C-^s - (V#e) /A] , (8) 

where C _1 is the error covariance matrix. Nevertheless, it is 
well known that this is an ill-posed problem, meaning that 
the solution of Equation JSJ is unstable with respect to the 
errors in the data, i.e. a small error in the initial data can 
result in much larger errors in the desired solution. 

The second information is related to the physics, and 
it allows us to ensure that the solution is physically accept- 
able. We use an approach known as regularization, which 
involves introducing additional a priori information in order 
to solve an ill-posed problem. This information is usually of 
the form of a penalty for complexity, such as restrictions for 
smoothness, by reducing the number of degrees of freedom 
of the problem. It includes various facts known about the 
solution a priori, e.g. the total cumulative mass in clusters 
must be positive, and monotonically increasing regardless 
the observational data. 

It is then necessary to introduce a means of balancing 
the regularization constraints with the accuracy to which 
the solution reproduces the data, which is parameterized by 
X 2 - This has been accomplished by using the Lagrangian 
multipliers to find a solution; in particular we minimize the 
function: 

L(e,\)= X 2 (e) + \C(e) (9) 
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Figure 1. Illustration of the onion-skin model adopted for the 
geometrical deprojection. We assume that the cluster is spheri- 
cally symmetric and it has a onion-like structure, with n concen- 
tric spherical shells and m rings or annuli (n > m). The matrix 
V] is a rectangular m X n matrix, whose entries represent the 
volume of the j-th shell contained inside the i-th annulus (with 
j J* i). The SZ boundary (the most external annulus) will be rep- 
resented by the mth ring of area A m and radius R m . The shells 
(m + 1, m + 1, .., n) of radius R > R m (long-dashed line) account 
for the edge effect. 



where C(e) is a function that is minimized when the so- 
lution best satisfies the regularity constraint, and A is the 
smoothing parameter. 

We adopt the Philips- To wmey regularization method 
l|Phillipslll96a lT womeviri963l ). Generally, this method con- 
sists of minimizing the sum of the squares of the k th order 
derivates of the solution around each data point. For our 
purposes, we define the "smoothness" constraint as the min- 
imum of the sum of the squares of the second derivatives, so 
that: 

n-l 

C(f) = X>^-i + 2b i * ^ + c ^-+i) 2 = e # H # e ( 10 ) 

3=2 

where the a, = 1, fe, = — 1, Cj = 1, and the matrix H is 

defined in lTwomevI (|l963h . 

The solution of Equation @ reads (see iBouchetll 199ot ) : 

e = (V^C-^V + A .ff) _1 #V#C- 1 # (s ■ A) (11) 

By varying A, we can therefore vary the degree to which 
the solution is dominated by consistency with the data (i.e. 
it is statistically acceptable) or with the regularizing con- 
straints (i.e. it is physically acceptable). The choice of A is 
therefore critical to obtaining a reliable solution. To choose 
the value of A we use the cross-validation, which offers an 
objective and robust way of choosing the smoothing param- 
eter. Cross-validation consists of predicting each data point 
s by finding a solution using all the data but that point. 



Table 1. Specifications of the CCAT SZ cluster observations. The 
stated frequency refers to the central frequency of the bandwidth; 
the noise level accounts for the confusion limit by dusty star- 
forming galaxies and refers to an integration time of 16 hours. 



v beam rms 

(GHz) (arcsec) (/^Jy/beam) 



93 


27.0 


20.0 


148 


16.8 


31.2 


221 


10.8 


25.0 


272 


9.0 


26.2 


347 


7.2 


38.8 


406 


6.0 


75.0 



The best value of A is the one that predicts the observations 
best, i.e. the one which minimizes the average quadratic er- 
rors between the da ta s and their predicted values or the 
reduced \ 2 (see e.g.. lBouchetl fl995. for further details). 

It is worth emphasizing the improvement of this method 
with respect to the standard deprojection technique in X- 
ray spectra (or surface brightness), where customarily the 
desired solution is inferred via inversion of Equation (0 un- 
der the assumption that n = m, neglecting the edge effect. 
In this case V would be an upper triangular matrix, whose 
inversion would allow us to determine e via Equation ©. 
While this approach is commonly employed in the deprojec- 
tion of X-ray physic al observables, such as the X-ray temper - 
ature or brightness (Kr iss et al.lll983l ; iMorandi et al.ll2007l ). 
for SZ data it would not be satisfactory, because the contri- 
bution of the shells with radius greater than the SZ bound- 
ary is non- negligible (since the SZ effect is linearly propor- 
tional to the electronic density, unlike the X-ray brightness), 
and the measurement errors a are not taken into account in 
the determination of the expectation value of e. Alterna- 
tively, one could extrapolate the profile of s via an analyti- 
cal function out to R n > R m and consider an n x n upper 
triangular matrix V in order to infer e via Equation © : nev- 
ertheless, in this case the solution could be systematically 
biased, depending on the assumed extrapolated profile and 
on the choice of outermost bin. Note that in our approach we 
corrected for the edge effect without assuming any a priori 
profile of s beyond the SZ boundary. 



4 MOCK CCAT SIMULATIONS 

In this work we analyze mock CCAT simulations of galaxy 
clusters. We start by describing the CCAT telescope in £14.11 
and the simulated cluster sample in £14.21 We then describe 
how we generate SZ maps from these simulations ( H4.3p . 
including the impact of relativistic electrons ( £|4.4|l and the 
primary CMB anisotropies ( £|4.5ll . 

4.1 CCAT SZ Observations 

Cerro Chajnantor Atacama Telescope (CCAT) is 25-meter 
submillimeter telescope that will be located at 5600 m alti- 
tude on Cerro Chajnantor in the Andes mountains of north- 
ern Chile. CCAT will combine high sensitivity, a wide field of 
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Table 2. Properties of the simulated cluster sample. We report 
-R500 1 the true global mass- weighted temperature T mw of the clus- 
ters, the minor-major and the intermediate-major axis ratios of 
the ICM and DM. The last column refers to the centroid-shift for 
the x, y and 2 projections, respectively. 



object 


-R500 






- m 


m - 


- '/2 


w 




(kpc) 


(keV) 


(ICM) 


(tot) 


(x/y/z) 


CL101 


1657 


7.44 


0.95- 


-0.98 


0.79- 


-0.88 


0.011/0.016/0.012 


CL102 


1397 


5.63 


0.92- 


-0.93 


0.85- 


-0.95 


0.023/0.023/0.009 


CL103 


1420 


4.84 


0.91- 


-0.99 


0.83- 


-0.92 


0.017/0.050/0.055 


CL104 


1394 


6.61 


0.90- 


-0.95 


0.79- 


-0.91 


0.002/0.005/0.004 


CL105 


1347 


5.67 


0.95- 


-0.96 


0.80- 


-0.89 


0.015/0.005/0.013 


CL106 


1202 


4.54 


0.72- 


-0.85 


0.67- 


-0.88 


0.039/0.053/0.037 


CL107 


1089 


3.61 


0.76- 


-0.80 


0.72- 


-0.75 


0.029/0.022/0.034 


CL3 


1016 


3.37 


0.80- 


-0.88 


0.73- 


-0.85 


0.007/0.010/0.007 



view, and a broad wavelength range to provide an unprece- 
dented capability for deep, large area multicolor submillime- 
ter surveys. In Table [1] we summarize the specifications of 
the proposed CCAT telescope, including frequency channels, 
angular resolution, and sensitivities for the integration time 
of 16 hours. 

4.2 Simulated Clusters Sample 

We analyzed a sample of high-resolution hyd rodynamical 
simula tions of galaxy clusters formation from dNagai et al.l 
l2007al lbl. he reafter N07a,b) which were performed u sing the 
ART code (jKravtsov et al.ll2002l : Iriudd et al.ll200ct ). In the 
present work, we analyze the z~0 outputs of the simulations 
that include radiative cooling, star formation, metal enrich- 
ment and stellar feedback at low redshift. In particular, we 
focused on a subsample of 8 clusters with the global mass- 
weighted temperature T mw > 3 keV, where the relativistic 
corrections would be large enough to robustly measure the 
cluster gas temperature from multifrequency SZ-only data. 
Each cluster is simulated using a 128 3 uniform grid with 8 
levels of refinement. Clusters are selected from 120/i _1 Mpc 
computational boxes, achieving peak spatial resolution of 
~ 3.6/i _1 kpc. The dark matter particle mass in the region 
surrounding the cluster is 9 x 10 8 h -1 Mq, while the rest of 
the simulation volume is followed with lower mass and spa- 
tial resolution. We refer readers to N07 for the details of 
these simulations. 

In order to assess the effect of cluster shapes on the 
derived 3D profiles, we calculated the intermediate-major 
and minor-major axis ratios of both the DM and the ICM. 
Assuming an onion-like density distribution in concentric 
ellipsoids, it is possible to det ermine the axial rat ios by di- 
agonalizing the inertia tensor (jSplinter et al.lll997l ): 

I', (i,j = 1,2,3), (12) 

where the sum is over all points xi = x, X2 = y, £3 = z of 
density p. With this definition we have the minor-major and 
intermediate- major axis ratios 7/1 and 772, respectively: 

Hfe)"* (13) 



with 771^*72^1, and where t/xx, U yy , and (7 Z z are 
the principal components of the diagonalized tensor, with 
(7 ZZ ^ Uyy ^ (7xx. We considered both the cases where 
p = picm and p = pt t , for the ICM and DM principal axis 
ratios, respectively. We also point out that some of the clus- 
ters show visible substructures, namely the objects CL103, 
CL106 and CL3. The most prominent substructures have 
been identified by visual inspection and masked out before 
the analysis. 

We then evaluated the SZ centroid-shift w of the Comp- 
ton parameter image, which is a proxy of the state of relax- 
ation and level of substructures of a cluster, w was deter- 
mined in a series of circular apertures centered on the cluster 
SZ peak out to the SZ boundary radius, roughly correspond- 
ing to 7?5oJE|. u> was defined as the standard deviation of the 
projected separations between the peak and centroid in units 
of 7?5oo- Table [2] presents the list of objects analyzed in the 
present paper, including R500, their temperature, principal 
axis ratios of both ICM and DM, and the measured centroid 
shifts. 

4.3 SZ maps 

We generate mock CCAT observations of the simulated clus- 
ters at 6 bands between ~ 100-800 GHz and with the beam 
size and expected noise level for 16 hours of CCAT observa- 
tions. We used the outputs of the simulated galaxy clusters 
at z — and placed it at the observing redshift of z = 0.2. 
Specifically, we simulated the SZ intensity maps AI(v) at 
each frequency by means of Equation {TJ and by using the 
3D gas pressure, temperature and peculiar velocity for the 
simulated galaxy clusters. For each cluster, the mock CCAT 
data (i.e., the SZ brightness AI(v)) are created for three or- 
thogonal projections by projecting both thermal and kine- 
matic SZ signals along the x, y, and z coordinate 
counting for the line-of-sight velocity structures (v x , v y and 
v z ) in simulations. 

Figure [2] presents the SZ brightness maps of CL101 (re- 
covered by applying Equation [4]) at each CCAT band. The 
size of the region shown in the figure is ~ 5.7 Mpc with 
our cosmological parameters, with a pixel size of ~ 22 kpc, 
and the region is centered on the minimum of cluster po- 
tential. We also model measurements at a given resolution 
by convolving the maps with a Gaussian beam profile of ap- 
propriate width corresponding to the resolution of CCAT at 
each frequency. 

4.4 Relativistic Electrons 

The presence of significant energetic electron populations in 
many clusters has been known from measurements of dif- 
fuse radio emission and non-thermal X-ray emission from 
relativistic electrons in several nearby clusters. The range 
of electron energies implied from these measurements is 
roughly 1-100 GeV. This can have appreciable impact on 
the reconstruction of the physical properties of the ICM 
from the SZ data, especially on the measurements of the 
cluster's peculiar velocity discussed in § 15.21 

8 -R500 is the radius within which the enclosed mean density is 
500 times the critical density of the universe. 
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93 GHz „„ 148 GHz „„ 221 GHz 




Figure 2. The SZ brightness maps of the cluster CL101 viewed along the z projection at 6 CCAT bands. The size of the region 

shown in the figure is 5.7 Mpc centered on the minimum of cluster potential, with a pixel size of 22 kpc. 



To assess the magnitude of this effect, we have included 
in the simulations the additional intensity change in the 
CMB spectrum due to non-thermal electrons: 



AJ rc l(z/) 

Io 



(e* - 1) ; 



(14) 



where r rc i = err / ?ie,rei dz is the optical depth of the rela- 
tivistic electrons. 

For simplicity, the non-thermal emission region has been 
assumed to have spherical symmetry with an electron spec- 
trum: 



N e (p,r) = k p 31 ■ g e (r) p > pi . (15) 

Here we assume si = 3, pi = 10, fco = 3.7 x 10~ 3 cm -3 
a radial distribution of the electron population given by: 



and 



9e(r) = 



1+1 - 



(16) 



with r c — 10 kpc and g e = 0.5. With this choice of feo, the 
non-thermal electron energy constitutes <C 1% of the thermal 
energy in the cluster, in agreement with the observations of 
Coma cluster (|Shimon fc RephaelJl2002l '). 

Therefore, the contribution of the non-thermal SZ signal 
is generally small compared to the tSZ signal (A7 ro i/AJ t h <C 
10 -2 ) in most bands, except near the crossover frequency at 



218 GHz (e.g., A/ re i/A7 th < 20% at 211 GHz). Hence rela- 
tivistic electrons have a negligible impact on the cluster tem - 
perature determination (see also, IShimon fc Rephaelj|j)02h . 
but they can lead to significant biases on the peculiar veloc- 
ity determination, which stems mainly from SZ measure- 
ments around 218 GHz. Note that, even if we increase fco by 
a factor of five in order to mimic the p roperties of clusters 
hosti ng strong radio halos (e.g. A2199, Shimon fc Rephaelil 
2002), we still find a small impact of the relativistic electrons 
on the derived cluster temperature. 



4.5 Primary CMB anisotropies 

Primary CMB anisotropies are another potential source of 
contamination in SZ observations. CMB anisotropies, AT, 
translate into uncertainties on the SZ signal AI(v) by the 
derivative of the blackbody spectrum with respect to the 
temperature, i.e. AI(v)/I = x A e x /{e x - l) 2 AT/T cmb . To 
assess the effect of the primary CMB anisotropies, we first 
produce a realistic CMB map T cm b + AT in order to generate 
the AI(v)/Io simulated images described in i]4.3l We then 
explicitly includ in the covariance matrix C (see 331) both 
the uncertainties due to the CMB correlated signal and the 
thermal noise of the instrument as: 



C = N + C cn 



(17) 
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Figure 3. The temperature (upper panel) and pressure (lower panel) profiles retrieved from mock CCAT analyses of CL101 and CL104. 
The black-solid represents the true mass-weighted temperature (pressure) of the clusters in the ith spherical shell, whose width is 
graphically represented by the horizontal solid lines. The open points with errorbars are the temperature (pressure) as retrieved from 
the multifrequency SZ-only analysis for the z projection, while the triangles (squares) refer to the x (y) projection (errorbars are omitted 
for clarity). We also report the fractional total scatter at the bottom of each panel. 



where C cm b is the CMB covariance matrix, N is the noise 
diagonal matrix u~ 2 , and a is the thermal noise of the in- 
strumental. 

We simulated realistic maps of the CMB given an 
input power spectrum C acco rding to the WMAP 7- 
year results |jarosik et al. l201Cf ). We then evaluated the 
product of the WMAP power spectrum Ci and the 
beam transfer function ~ exp(— 0.5 1 (I + 1)# 2 ), with 
I = 7r/arctan[0.5-v/(M 2 + « 2 )] ~ V 27V (« 2 +« 2 ), 9 = 
/ \/8 log (2), and 0f w hm is the beam of the telescope. 

Given a CMB image AT(x,y) characterized by N x N 
pixels and angular resolution Ax = Ay, we performed 
Monte-Carlo simulations of AT(x, y) in the Fourier domain 
by creating a complex visibility V(u, v) with Ait = Av — 
1/(N Ax). For each point in the Fourier domain we com- 
pute its distance I and the corresponding C'i from the input 
power spectrum. We set the expectation value and variance 
of V(v, u) to be zero and \fClAu, respectively. We then per- 
form a Gaussian randomization of both the real and imagi- 
nary parts of V(u, v), taking into account its conjugate sym- 
metry. We set the central visibility V(0, 0) to be zero, so the 
sum of the CMB fluctuations AT is zero over the whole im- 
age. We then obtain realization of the CMB map AT(x, y) 
by performing the inverse Fourier transform of V(u, v). 

In order to assess the impact of primary CMB 
anisotropies, we employed a lead— trail (L— T) differencing 
scheme where two fields (the source and the CMB fields) 
are observed in succession at the same hour angle and sub- 
tracted from each other. This removes the ground signal 
and any other potential spurious signals. The angular dis- 
tance between the L— T fields has been assumed to be 45 
arcmin. The impact of the primary CMB anisotropies de- 



pends mainly on the angular resolution (beam), the angular 
size of the cluster, the width of the annuli where the signal 
ty(i>; T(r), /3(r)) is averaged out and deprojected, and the 
angular distance between the L— T fields. The contribution 
of the primary CMB anisotropies increases toward larger an- 
gular scale out to about one degree, and it becomes a domi- 
nant source of uncertainties beyond the SZ boundary of the 
observation on scale > 4 — 6 arcmin, roughly corresponding 
to 7?5oo- This also sets the outer radius out to which the 
three dimensional profiles can be reliably measured using 
the multifrequency SZ observations. On scales < 1 arcmin 
the CMB anisotropies is negligible compared to the thermal 
noise. 



5 RECONSTRUCTION OF 3D CLUSTER 
PROFILES 

5.1 Gas temperature, pressure and total mass 

We test the performance of our inverse method for recon- 
structing the desired 3D properties of the ICM using hy- 
drodynamical simulations of galaxy clusters. Using the non- 
parametric method, we deproject s in order to infer e (33J) , 
by rebinning the SZ brightness in circular annuli around 
the cluster center for each frequency and projection axis. 
Given the spectral signature of 9(u; T(r), /3(r)), we then sep- 
arate the physical variables P(r),T(r),/3(r) by comparing 
the deprojected observed quantity e to its theoretical pre- 
diction. Errors have been evaluated via Monte-Carlo ran- 
domization, by generating correlated random numbers p 
with probability density function / = /(s, C) where s is 
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Figure 4. The recovery of 3D mass profile from mock CCAT analyses of CL101 and CL104. The solid line represents the hydrostatic 
mass using the true mass-weighted temperature and pressure, while the dashed one is the true spherically-averaged mass. The open 
points with errorbars are the recovered hydrostatic mass along the z projection, while the triangles (squares) refer to the x (y) projection 
(their errorbars are omitted for clarity). 



the expectation value and C is the covariance. The final er- 
ror estimate is accomplished by diagonalizing the covari- 
ance matrix C. In practice, by using its eigenvectors V , we 
rotate into a space in which the covariance matrix becomes 
C' = diag((Ji , erf , ), i.e. it is diagonal with eigenvalues 

erf, erf') •••! Op'- We can thus treat the rotated parameters p' 
as independent and uncorrelated measurements with proba- 
bility density / = f(s, C'). To obtain the probability density 
function / = f(s, C) we use the following relation: 

P = V#p' . (18) 

Figure [3] shows the recovery of 3D temperature and 
pressure profiles from the mock CCAT analysis of CL101 
and CL104 clusters viewed along three orthogonal projec- 
tions. CL101 is a massive, dynamically active cluster, which 
experiences violent mergers at z ~ 0.1 and z ~ 0.25. CL104 
is a similarly massive cluster, but with a more quiescent his- 
tory. This cluster has not experienced a significant merger 
for the past 6 Gyrs, making it one of the most relaxed sys- 
tems in the N07 sample. 

We find a remarkable agreement between the "mea- 
sured" temperature (pressure) profiles and the true mass- 
weighted temperature (pressure) for the relaxed cluster 
CL104. We find a variation of order 5 — 15% in both the 
recovered pressure and temperature profiles obtained for dif- 
ferent projections of the same cluster. The recovery of the 
ICM profiles is still quite good for the unrelaxed cluster 
CL101 with the systematic bias at the level of 15 — 30%, 
due to the combination of the presence of substructures and 
the departure of the ICM distribution from the spherical 
symmetry. Moreover, internal bulk motions can introduce 
additional biases in the derived ICM profiles, especially in 
disturbed systems (e.g., characterized by the high values of 
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Figure 5. Upper panel: ratio between the measured and true 
mass-weighted temperature averaged inside a sphere of -R500 for 
the sample of simulated clusters. Middle panel: ratio between the 
measured and true Compton parameter integrated within -Rsoo- 
Lower panel: ratio between the measured and true peculiar veloc- 
ity along the line-of-sight, mass- weighted averaged inside a sphere 
of radius of Rsoo- For all the panels, the open points with error- 
bars refer to the z projection, while the triangles (squares) refer 
to the x (y) projection (errorbars are omitted for clarity). 
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Figure 6. Reconstruction of 3D temperature profiles from the mock CCAT analyses of CL101 and CL104 viewed along the z projection, 
in which we add contaminants to the tSZ signal one at a time, starting with the kSZ signal, the relativistic electrons, the primary CMB 
anisotropics, and the thermal noise in order. The different colored points reflect different astronomical and instrumental effects which 
degrade the recovery of profiles (see legend in the figures). The black-solid line represents the true mass- weighted temperature of the 
clusters in the ith spherical shell, whose width is graphically represented by the horizontal solid lines. We also report the fractional total 
scatter at the bottom of each panel. 



w in Table [2j. Note that via our non-parametric method 
the full information about the ICM profile is preserved in 
converting the SZ brightness profile to 3D temperature and 
pressure profiles. 

We also assess the impact of the primary CMB 
anisotropics on the derived ICM profiles, by comparing the 
recovered profiles from the maps with and without CMB 
anisotropies, while always including the thermal noise. We 
find that the bias due to the primary CMB anisotropies is 
small in the cluster cores, but they become comparable to 
that of the thermal noise in the outer radii (r > i?5oo). 

In Figure [4] we show the mass profiles obtained by ap- 
plying the hydrostatic equilibrium equation on the temper- 
ature and pressure profiles derived from the mock CCAT 
analysis. Similarly, we also compute the hydrostatic mass 
profiles using the mass-weighted temperature and pressure 
measured directly in simulations. We then compare these hy- 
drostatic mass profiles to the true mass measured directly 
from simulations in order to assess the bias in the recov- 
ered mass profiles. We find that the multifrequency CCAT 
SZ analysis of individual cluster can recover the hydrostatic 
mass profiles very accurately in the entire radial range for 
the relaxed cluster. Note, however, that the hydrostatic mass 
is biased low compared to the true mass by 5 — 10% (N07b), 
primarily due to the violation of the hydrostatic condition 
even in the inner regions of a very relaxed cluster (CL104). 
In the unrelaxed cluster (CL101), the bias is larger espe- 
cially in the outer regions where the departure from the 
hydrostatic equilibrium is even more significant. Note that 
the bias in the hydrostatic mass could be corrected if we 



could measure the level of bulk and turbulent gas motions 
in clusters (e.g., by kSZ effect) and account fo r the non- 
thermal pressure provid ed by these components l|Lau et al.l 
120091 ; iNelson et al]|2012r i . 

In Figure [5] we demonstrate the recovery of the global 
properties of the ICM, by showing the ratios of the measured 
and true mass-weighted temperature, integrated Compton 
parameter (Ysz), mass- weighted line-of-sight peculiar veloc- 
ity within a sphere of radius R^oo for the sample of simulated 
clusters. We find a non- negligible scatter in both the recon- 
structed temperature and Ysz- For 16 hours of CCAT obser- 
vations, the total scatters (intrinsic plus statistical measure- 
ment errors) are ~13% for the mass-weighted temperature 
and 4% for Ysz, while the intrinsic scatters (variance of the 
data points without the measurement errors) are relatively 
small (~ 4% and ~ 3%, respectively), primarily due to sub- 
structures and triaxiality of the ICM distribution. 

In Figure [6] we assess the impact of various sources of 
uncertainties by adding contaminants to the tSZ signal one 
at a time, starting with the kSZ signal, the relativistic elec- 
trons, the primary CMB anisotropies, and the thermal noise 
in order. We find that the temperature profile is recovered 
very well for CL104. The recovery is not as accurate for 
CL101, in which the recovered temperature could be biased 
by 5 — 25% at any given radius. In absence of other sources of 
noise (i.e., using the tSZ data alone), the reconstructed 3D 
temperature profiles are biased by only ~ 4% for CL104 and 
~ 9% for CL101. Biases are generally larger for the unre- 
laxed systems that exhibit more significant departures from 
the spherical symmetry indicated in Table [2] At large radii 
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Figure 7. Absolute value of the ratio between the kinetic and 
thermal SZ effect via an observation of CL101 with a receiver 
of bandwidth equal to 7.5 GHz and centered on 218 GHz. We 
take a Hanning window for the response of the radio receiver and 
we projected the datacube along the z direction. The size of the 
region shown in the figure is about 5.7 Mpc with our cosmological 
parameters, with a pixel size of about 22 kpc, and the region is 
centered on the minimum of cluster potential. 



(r J> R500), we find that the primary CMB anisotropies are 
one of the main sources of systematic uncertainties in the 
3D temperature reconstruction. For example, at r = 3i?soo 
the bias due to the parimary CMB anisotropies is ~ 10% of 
the SZ signal at 93 GHz. In the inner regions (r < -R500), 
the primary CMB anisotropies introduce the scatter at the 
level of 7% and 10% for CL104 and CL101, respectively. Our 
analyses indicate that the relativistic electrons, peculiar ve- 
locity and thermal noise have little effect on the recovered 
3D temperature of a hot (T mw .> 3 keV) galaxy clusters. 



5.2 Cluster peculiar velocities 

Galaxy clusters are excellent tracers of the large-scale veloc- 
ity field in the universe. The kSZ effect, due to the relative 
motion of the cluster with respect to the rest frame of the 
CMB, enables direct measurements of the cluster peculiar 
velocity independently of the cluster distance. In principle, 
this technique can be used to measure deviations of the clus- 
ter's motion from the Hubble flow with high precision at 
any redshift. This is in contrast to the standard candle (e.g., 
SNe) whose measurement errors increase monotonically with 
redshift. In this section we evaluate how well we can mea- 
sure the peculiar velocity of clusters using multifrequency 
SZ observations. 

In practice, there are a number of complications for 
the kSZ measurements of the cluster peculiar velocity: (i) 
the kSZ effect is much smaller than the tSZ effect and is 



entangled with other contributions to Comptonization of 
the CMB, such as energetic non-thermal electrons: difficul- 
ties might arise in separating these different spectral signa- 
tures in SZ signal P(r) ^{v; T(r), /3(r)) via multifrequency 
SZ data; (ii) the blackbody kSZ distortion is degenerate with 
the spectrum the primordial CMB anisotropies, which is 
also a blackbody; (iii) the observed kinematic SZ signal gets 
significant contributions from the internal gas motions in 
galaxy clusters, limiting the accuracy of the peculiar veloc- 
ity meas urement at the lev el of 50 — 100 km/s for individual 
clusters jNagai et alj|2003f ). 

Concerning the first point, the cluster peculiar veloc- 
ity is often estimated from the kSZ maps alone, by focus- 
ing on observations at the crossover frequency at 218 GHz 
where the tSZ emission vanishes. However, the exact value 
of the crossover frequency depends wea kly on the electro n 
spectrum (gas temperature and density. lBirkinshawlll999l ). 
Moreover, real observations are obtained with receivers with 
a finite bandwidth typically of order a few GHz. The kSZ 
observations are contaminated by other signals, including 
tSZ emission and non-thermal SZ from relativistic electrons. 
This is illustrated in Figure[7] which shows the ratio between 
the kinetic and thermal SZ effect obtained by receivers with 
bandwidth of 7.5 GHz and centered on 218 GHz. We can 
see that the thermal and kinetic SZ signals are comparable 
in intensity; i.e. the tSZ effect does not fully vanish because 
of the finite width of the bandwidth and the dependence of 
the crossover frequency on the gas temperature and density. 

To address these issues, we analyze mock CCAT simula- 
tions in order to assess how well we can measure the peculiar 
velocity of galaxy clusters. Figure [S] shows the reconstructed 
streaming velocity profile /9(r) as a function of the cluster- 
centric radius. We find that the velocity radial profile /?(r) in 
the ICM matches the true spherically-averaged cluster pe- 
culiar velocity. Our retrieved streaming velocity profiles sug- 
gest that systematic uncertainties in the measured /3(r) is 
typically of order of 100 km/s, but there are occasional out- 
liers with bias of order a factor 2 in a certain projection axis 
of disturbed clusters. For example, the largest bias is found 
in CL106, which is a disturbed system with a large value of 
w (Table [2}. Therefore, one way to control such a large bias 
is to consider only relaxed systems with a low value of w 
and average the peculiar velocity out to the SZ boundary. 
Note also that the biases are present even in relaxed systems 
due to frequent minor m ergers that continuously perturb th e 
clusters gas atmosphere (|Vazza et a l. 2009; Lau et alT 2009). 

Next, we consider the mass-weighted peculiar velocity 
averaged within a circular aperture. Averaging out within 
a fixed radius has two advantages: (i) the mass-weighted 
peculiar velocity is more easily obtainable with SZ obser- 
vations with moderate angular resolution; (ii) the bias in 
the measured mass-weighted streaming velocity tends to be 
smaller than the differential radial profiles /3(r). In Figure [5] 
we present the ratios between the measured and true pecu- 
liar velocities obtained along three orthogonal projections, 
where the recovered velocity is mass-weighted averaged in- 
side a circular aperture of radius equal to the boundary of 
the observation. For most of the objects, the average pe- 
culiar velocities derived from different projection axes are 
consistent within statistical errors of order 10 — 100 km/s 
(corresponding to a fractional uncertainty of ~ 10 — 15%). 
This dispersion sets a theoretical lower limit on the accu- 
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Figure 8. The recovered spherically-averaged streaming velocity profiles /3(r) from mock CCAT analyses of simulated clusters viewed 
along three orthogonal projections (x, y, z). The black-solid lines represent the true spherically averaged mass- weighted streaming velocity 
of the clusters in the ith spherical shell, whose width is graphically represented by the horizontal solid lines. The points with errorbars are 
the streaming velocities retrieved from the mock analysis. The dashed lines refer to the true mass-weighted streaming velocity averaged 
inside a circular aperture of radius equal to the boundary of the observation. 



racy with which cluster peculiar velocity can be measured 
for an individual cluster. The bias in the averaged peculiar 
velocity mainly arises from: (i) internal motions of the ICM, 
which gives an sy stematic error of around 50 — 100 km/s 
l|Nagai et al.ll2003h ; (ii) triaxiality, which introduces an sys- 
tematic error of <J 50 km/s; (iii) systematic intensity change 
A/ ro i(^) due to relativistic electrons, which gives an sys- 
tematic error of < 20 km/s. These errors are also largely 
independent of the cluster redshift, in contrast to galaxy pe- 
culiar velocity surveys. Therefore, kSZ measurements of the 
cluster peculiar velocities might serve as promising tools for 
studying cosmology, including the expansion history of the 
universe, the matter density of the universe, the primordial 
power spectr um normalization, and the dark energy equa- 
tion of state (jBhattacharva fc Kosowskvll2007h . 



6 DISCUSSIONS AND CONCLUSIONS 

In this paper we presented a non-parametric inversion 
method for reconstructing 3D temperature and pressure pro- 
files of the ICM from multifrequency SZ observations. In 
order to assess the accuracy to which the 3D physical pa- 
rameters can be recovered from observations, we analyzed 
mock Cerro Chajnantor Atacama Telescope (CCAT) obser- 
vations of simulated clusters, while taking into account sev- 
eral sources of uncertainties associated with instrumental 
effects and astrophysical foregrounds. 

We showed that our method enables the recovery of 
the 3D ICM temperature and pressure profiles out to large 
radius (r ss R500) of a hot (Tx J> 3 keV) cluster, while 



retaining full information about the gas distribution. Ac- 
curate recovery of the ICM pressure and temperature pro- 
files, in turn, enables accurate determination of the mass- 
weighted temperature, integrated SZ signal (Tsz), and the 
cluster mass profiles assuming the hydrostatic equilibrium 
of the gas in the cluster potential. We found that the bias in 
the reconstructed 3D ICM and mass profiles are dominated 
by substructures and asphericity of the cluster gas distribu- 
tion, while the contaminations by thermal noise, relativistic 
electrons and the CMB anisotropics are small. Note also 
that the hydrostatic mass is a biased estimate of the true 
cluster mass, due to the presence of non-thermal pressure 
provided by bulk and turbulent gas motions. A joint analy- 
sis of SZ and lensing measurements may help constrain the 
biase s due to triaxial shape and non- thermal pressure fur- 
ther l|Morandi et alj|201(l l201ll . |2012h . 

We also assessed the accuracy to which the cluster's 
peculiar velocity could be measured with upcoming multi- 
frequency observations of the kSZ effect. We find that there 
is a non-negligible total scatter of ~ 10 — 15% for individual 
clusters. For 16k hours of CCAT observation, the measure- 
ment errors are still dominated by statistical uncertainties, 
instead of systematic uncertainties. The systematic bias in 
the averaged peculiar velocity mainly arises from: (i) inter- 
nal motions of the ICM (w 50 — 100 km/s); (ii) triaxiality 
(<C 50 km/s); (iii) systematic intensity change AI Te i(v) due 
to relativistic electrons (< 20 km/s). Thus, CCAT measure- 
ments of the cluster peculiar velocities via kSZ effect should 
serve as unique and potentially powerful probes of the dy- 
namics of the universe and cosmological parameters. 

The inversion method presented here should be a valu- 
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able tool for analyzing multifrequency SZ observations of 
galaxy clusters. While the SZ effect has been used mainly for 
constraining the electronic pressure to date, high-resolution, 
multifrequency SZ observations (e.g., by CCAT) will enable 
detailed characterization of the temperature and velocity 
structures in the ICM, without follow-up X-ray observa- 
tions. SZ observations are especially unique and powerful 
in probing the outskirts of high-redshift clusters and com- 
plementary to X-ray observations which are sensitive to the 
inner regions of nearby clusters. In this sense, CCAT SZ 
observations are complimentary to the upcoming eROSITA 
X-ray cluster surve>|j with comparable angular resolutions. 
Furthermore, CCAT SZ observations will provide unique 
constraints on constraints on the cluster peculiar velocities 
and possibly internal gas flows especially in the outskirts of 
high-redshift clusters, which are difficult to measure using 
the upcoming Astro-H missio via shifting and broaden- 
ing of X-ray spectral lines. The present work is a step for- 
ward toward robust reconstruction of the 3D ICM, mass, 
and streaming velocity profiles for a sample of clusters and 
promises to help advance cosmological studies based on the 
cluster growth over cosmic time in the coming decade. 
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